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Abstract 

The dispersion relations and ^-matrix of the one-dimensional Hubbard model at half 
filling are considered in a certain scaling limit. (In the process we derive a useful small- 
coupling expansion of the exact lattice dispersion relations.) The resulting scattering 
theory is consistently identified as that of the SU(2) chiral-invariant Thirring (or Gross- 
Neveu) field theory, containing both massive and massless sectors. 



1. Introduction 

The Hubbard model [1] describes electrons on a lattice with on-site interaction only, 
in addition to a standard nearest-neighbor hopping term. In two dimensions the model 
has received much attention lately in connection with high-T c superconductivity. Some 
of its properties are believed to be similar to those exhibited by the one-dimensional 
model, which is exactly solvable by means of the Bethe Ansatz technique [[| . In this paper 
we discuss certain aspects of the scaling limit of the one-dimensional model, which are 
relevant for its large-distance asymptotic behavior. 

The hamiltonian of the linear Hubbard model is given by 

^ = "2EE (4,c i+v + C ; +Ii/J> ) - 2[/^(n i , T - ^(n^ - -) , (1.1) 

j=l <r=T,4 j = l 

where c J)(T are canonical fermionic annihilation operators, j labels the sites of a periodic 
chain of length L (which is taken to be even), a labels the two spin degrees of freedom, 
U is a real coupling constant, and rij i(7 = cj a Cj )(7 is the number operator for spin a on 
site j. (The overall normalization of H chosen in ( |1 . 1| ) will be convenient later on.) Since 
H commutes with the total number operator Y^j=i X/ CT =t l n i> cr ' ^ can ^ e diagonalized 
separately in eigenspaces of fixed number of "electrons" N. 

The model hasS an SO (4) = SU(2) x SU(2)/Z 2 symmetry, the two 5*7(2) 's pertaining 
to spin (s) and charge (c). The spectrum is built out of four fundamental excitations (alias 
quasiparticles) , forming a "spinon-antispinon" 5'L r (2) s -doublet and a "holon-antiholon" 
S'?7(2) c -doublet. However, this separation of spin and charge seen in the quantum numbers 
of the quasiparticles does not mean that the theory decouples into a tensor product of 
two S'?7(2)-symmetric models. For instance, there is a selection rule which allows only 
representations with integer total SU(2) S and SU(2) C spin in the spectrum (implementing 
the Z2 quotient which reduces the symmetry from the naive SU(2) x SU{2) down to 
SO(4)).@ 

The situation at half-filling N = L is of special interest. In this case one of the two 
quasiparticle doublets develops a mass gap while the other remains massless. The fact 
that the mass gap vanishes as the coupling tends to zero opens up the possibility for the 

1 For details on the following features of the model see [Q and references therein. 

2 In particular, this restriction implies that there are no single-particle states in the spectrum 
and thus the fundamental quasiparticles are "confined" . 
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existence of a scaling limit in which both massive and massless excitations survive in the 
spectrum of the resulting field theory. Our aim is to explore this possibility, which we 
will do at the level of the dispersion relations of the quasiparticles and their scattering 
amplitudes. The S'-matrix theory obtained this way is then identified as that of the SU(2) 
chiral-invariant Thirring (or Gross-Neveu) model, whose lagrangian is given by [5-7] 

C = i#pif> + g[$il>) 2 ~ (^7W] 

1 3 (12) 

a=l 

where if) is a doublet of Dirac spinors, 7^ are Dirac matrices in 1+1 dimensions, and a a 
are the Pauli matrices (the equality between the two lines of (|1.2 ) can be established with 
the aid of identities listed in the appendix of ||). 



The field theory (1.2) and its ^-matrix have been discussed in [8-16]. It is known, 
using bosonization, that the theory essentially decouples into a massless and a massive 



sector. (This statement holds modulo certain orbifolding, cf. |I5] ||17|1 , which is reflected for 
instance by "kinky" restrictions |18[ on the multiparticle spectrum.) The massless sector 
is described by the level one SU(2) WZW conformal field theory. The massive sector, 
on the other hand, can be viewed as a marginally relevant 5'L r (2)-preserving (integrable) 
perturbation of another copy of the same conformal field theory, where mass is generated 
dynamically through "dimensional transmutation" (for a construction of this sector of the 
theory from a scaling limit of the XXZ spin chain cf. [|l9| ).i 

The emergence of the SU(2) chiral-invariant Thirring field theory from the scaling 
limit of the half-filled Hubbard model, as described in the sequel, is not surprising. It 



was already noted on the basis of renormalization group and symmetry arguments in [15[] , 
where the continuum (low-energy) limit of ( |1 . 1| ) was considered. Nevertheless, we think 
that our complementary analysis is worthwhile. Related work can be found in ||20|| . 

The rest of the paper is organized as follows. In sect. 2 we define the scaling limit 
and derive the scaled dispersion relations. In the process the familiar Hubbard model 
dispersion relations (2.1)-(2.2) are rewritten in the form (2.3)-(2.4), which is most useful 
for analyzing the rather singular zero-coupling limit; we find this apparently new form, 
whose derivation is presented in the appendix, interesting by itself. The scaled S'-matrix 
obtained in sect. 3 is discussed in the final section. 



3 Somewhat confusingly, in the literature the S-matrix of the massive sector alone is occasion- 
ally referred to as that of the full theory (|1.2|). 
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2. The scaling limit 

We restrict attention to the attractive regime U > 0, where the spin excitations 
(spinons) are massive while the charge excitations (holons) are massless. The repulsive 
regime is dual to the attractive one in the sense that the properties of the two excitations 
are interchanged || [^TJ . The spin- wave dispersion relation is given in the parametric form 
(see 0@ and references therein) 



p j k -\ = k _ I ' du Jo(u)sm(u8mk) c _ wU 



n OO 

e s (k) = U — cosk + 

Jo 



cko Ji(u) cos(uj sin k) _ LoU 



(2.1) 



uj cosh.(uiU) 

where the J u {oj) are Bessel functions. The charge- wave dispersion relation reads 

f°° du J (a;)sin(u;A) f°° du J^uj) cos(u\) 

Pc(A) = - / , n , e c (A) = / — , (2.2) 

J uj cosh(o;f/) J uj cosh(ujU) 

where, for nonzero U, A G (—00,00). (In the free case U =0 one has |A| < 1, and, using 
formulas 6.693(1-2) of |22H , eq. (|2.2| ) reduces to p c (\)\ U Q = — arcsinA and e c (X)\ U = 
— A 2 , so that e c (p)\ UQ = cos p.) 

In order to take the scaling limit we need a more convenient form for the dispersion 

1\ 7T 

2> U 

as 



(2.3) 



relations. Let U n = (n+ Then, as shown in the appendix, eq. (|2.1|) can be expanded 

2 °° 1 

p s {k) = jf^2 jr K °^ Ur ^ sin H U n sink) 

U n=0 Un 

2 °° 1 

e s (k) = —y j —K 1 (U n ) cosh(?7 n sin fc) , 

n=U 

for k G (-■§■, §), while for |A| > 1 the following expansion of eq. ( |2.2| ) is valid: 
Pc(A) = sgn A (-l + l£^I (U n ) e -W\ 

n=0 



n=0 n 

Here K u (z) and I v (z) are the modified Bessel functions. 



(2.4) 



4 Our conventions differ from those in [fl by an overall factor of | in the hamiltonian (1.1) 



and a change in the sign of U. 
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From eq. ( |2.3| ) we read off the spinon mass gap (cf. P[||23||) 



m) = e . {0) = i± , ,, 5) 

7T An + 1 



n=0 



which vanishes for smaii U according to pi 



A{U) ~ l^ e -w as t/^o+. (2.6) 

7T 

(For large U, on the other hand, the leading behavior is A(U ^> 1) ~ ^ Xl^Lo 2n+i = T' 
in between A(i7) increases monotonically for all U > 0.) The vanishing of the mass gap 
as U — > is a necessary condition for the existence in this limit of a nontrivial field 
theory where the spinon survives as a massive particle of finite mass M. Introducing 
a dimensionful lattice spacing a (so that H of is replaced by H/a), we define the 
corresponding scaling limit by 

a,U^0 + with M= lim = lim l^ e "w fixed. (2.7) 

a, {/— »o a a,C7— >-o 7ra 

To obtain a whole massive dispersion curve in the scaling limit, the lattice rapidity 
variable k has to be rescaled according to 

k -> such that = is finite. (2.8) 



As a result, it follows that only the n=0 term in ( |2.3| ) survives in the limit ( |2.7| ), leading 
to the scaled momentum and energy 

P s = lim^^ = Msinhfl , E s = lim = Mcoshfl . (2.9) 

a a 



This is a relativistic massive dispersion relation E S (P) = \J M 2 + P 2 in 1+1 dimensions, 
parameterized in terms of the customary rapidity variable 9 in the continuum. Note that 
the overall normalization of the hamiltonian ( |1 . 1|) is such that the "speed of light" is 1. 
Turning to the ho Ion dispersion relation, we see from (|2.4| ) that the energy e c (A) van- 



ishes as A — > ±oo (for fixed U), the leading behavior being e c (|A| 1) ~ 2Ii(-^j)e~ 7r \ x ^ 2U . 
The momentum p c (A) approaches in this limit, like p c (|A| 3> 1) ~ sgn A(— ^ + 
2/o(^/)e _7r ' A '' /2 ' 7 ). This implies that the dispersion curve linearizes around these points: 

^(p-T?) ~ v c {U)C--\p\) , v c (U) = . (2.10) 



The speed of the charge wave v c (U) was first given in [E3 (cf. also [25]). We note that 
v c (U — > + ) = 1, already suggesting that in the scaling limit the ho Ions become massless 
particles, traveling with the same "speed of light" as obtained from the massive dispersion 
relation of the spinons. 

To obtain the full massless dispersion relation in the scaling limit Q2.7T) , we first note 
that as U — > the holon dispersion relation (|2.4j) becomes linear for all |A| > 1, and not 



just as |A| — > oo. In fact, we may define the rescaled rapidity variable f3 by letting 

A -> ±2 such that (3 = ± ^—^j-^ is finite, (2.11) 



when the limit (|2.7| ) is taken. Here and below the upper and lower sign choices apply to 
right- and left-moving excitations, respectively, which have to be treated separately. Now 
using the asymptotics of the functions I„(z), one finds from ( |2.4j) and ( |2.11| ) that the scaled 
dispersion relation takes the form 

P c = Kmi(»,(A)±f) = ±fe« , E c = lim^ = ^e±" , (2.12) 

in the scaling limit ( |2.7| ). 

Eq. (|2.12 ) provides the standard parameterization [26-28] of the dispersion relation 



E C (P) = \P\ of a massless particle in a (1+1) -dimensional quantum field theory, with 
(3 G (— oo, oo) (for both right- and left-movers) and M being some mass scale. This mass 
scale is arbitrary and irrelevant if the theory is conformal, since no observables depend 
on it and a change in it can be absorbed by a redefinition of (3. We find it satisfying 
that in our case, where the theory does have a massive sector, this mass scale turns out 
to be exactly equal to the mass M of the massive excitation Q2.9Q , when one uses the 
most natural definition ( |2.11|) of (3. (Of course one can define (3 as in (|2Tfl) but with 



A — > ±(2 + const • U), which would have shifted (3 and hence rescaled M of ( |2.12j ) by a 
finite factor; the choice const=0 is what we call the most natural one.) 



3. The S'-matrix in the scaling limit 

In the previous section we found that the spectrum of the model in the scaling limit 
( |2.71) consists of two doublets (labeled by s and c) of particles, one massive and the other 
massless. We can now use the known S'-matrix computed in |4|] for the spin chain to obtain 
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the ^-matrix of the continuum theory, simply by reexpressing the amplitudes of [|J in terms 
of the rescaled variables (|2.8|), ( |2.11| ). 

Due to the SU(2) S x SU{2) C symmetry of the model, the two-particle S'-matrix is 
block-diagonal with four 4x4 blocks S xy (x,y G {s,c}), corresponding to the scattering 
sectors s-s, s-c, c-s, and c-c. (In each block the rows and columns correspond to incoming 
and outgoing particles; they are labeled by a pair of SU (2) quantum numbers, namely spin 
up '+' or down '— '.) In the spin chain, the amplitudes in each sector are functions of a 
single variable [i which is defined as 



sin ki — sin k 2 \ 

W 
sin k — A I 



p = — in s-s 



2U 
Ai — A2I 



in s—c or c—s) 



m c—c 



(3.1) 



2U 

where the index .7=1,2 refers to the two scattering particles. In terms of the rescaled 
rapidities (|2.8|) and (|2.11| ), where U — > + , kj — > 0, and Xj — > this becomes 

\6i — 6 2 \ 

[i — > m s—s 

7T 

^00 in s—c (3.2) 

-\Sl=M f or R-R or L-L :„ _ _ 

7r m c — c , 

00 for R-L or L-R 

where R and L indicate right- and left-movers, respectively, in the massless sector. 

Using the results of || we now find the following S'-matrix amplitudes in the scaling 
limit: 

s 8 s(di,e 2 ) = So(e)(-^—i--^-u) , (3.3) 

\t) — in t) — in J 
/ 14. i e nv \ 

S 8C {0,(3) = lim (-1 I . - - I = il, (3.4) 

^ R )(/3 l5 /3 2 ) = S^iPufo) = -S SS {(3 U (3 2 ) 

SS^Hfafo) = lim S (-n)I = il , (3 ' 5) 

where 6 = 61—62, I and IT are the identity and permutation matrices {i.e. I^ b b = 
8%' 8% and Ii a 2' = 8%8%' with a, 6, a', 6' G {+, -}), and 

f +f ri-^ P l Jo u coshc^ r 



(Of course Jq(0) = 1 on the most rhs of (3J5); it is inserted there, as well as the arbitrary 
variable u > which can simply be rescaled away, in order to exhibit an amusing relation 
between the phase shift asso ciated with S (9) and p a (k) of Q at k such that = 9, 



2 it 



4. Discussion 



Eq. ( |3.3| ) is identified as the ^-matrix of the massive sector of the 5'C/(2)-Thirring 
model [9-10,12-14]. It is equal to the limit g — > (— §) + of the ^-matrix of the ordinary 
massive Thirring model |29]]|3(]] (which, up to a sign ||17|| , is also that of the sine-Gordon 



model in the limit (3 — > v87r ). Eq. (|3.5|) , on the other hand, has been recently pro- 
posed |27j (cf. also JST|) to describe the massless scattering theory associated with the 
level one SU{2) WZW conformal field theory. Finally, the fact that S sc (9,f3) turned out 
to be rapidity-independent indicates that the massive and massless sectors essentially de- 
couple in the scaling limit. These observations are all in concert with the identification of 
the full scaled model as the SU(2) chiral-invariant Thirring field theory. 

Furthermore, as mentioned in the introduction, this field theory is not quite a true 
tensor product of two sub-theories, due to nontrivial "gluing" of sectors in their spectra. 
In the spirit of []I7fl, it appears that the phase factor i in (|3.4| ) signals this effect at the S- 



matrix level. It would be interesting to investigate what consequences this factor may have 
on the correlation functions of the theory.! However, for doing that a better understanding 
of scattering theories involving massless particles is needed, as well as an extension of the 
form factor bootstrap program to their framework. 



Acknowledgements. I would like to thank F. Efiler for discussions. This work was 
supported in part by the US-Israel Binational Science Foundation. 



5 Recall the work of [ {32f , where the spin correlation functions in the Ising field theory [33] were 
reconstructed using the form factor bootstrap program [14]. In this model, the nontrivial sign of 
the S-matrix S(0) = — 1 leads to correlators which are expressed in terms of solutions to Painleve 
equations |33[ | , rather than simple Bessel functions which arise in fermion correlators in the trivial 
theory of a free massive Majorana fermion whose S-matrix is simply S(6) = 1. Cf. also [34]. 
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Appendix A. Analysis of the dispersion relations 

We start with the massless case which is relatively simple. To derive eq. (|2.4|) from 
fl2.2p we first expand the Bessel functions J v (u) in powers of u and integrate the resulting 



series term by term, using formulas 4.111(3,4,7) of [22|. This gives 



Pc(A) = -sgn A 



7T 



+ 2 arctan e 2 



1 



2 A' 



2k-l 



1 



' 2 2fc (fc!) 2 V2E/7 ^dx / cosh 

2/s 



x 



(A.1) 



1 



fc=0 



2 2/ -+ :1 A-!(A- + 1)! i.2u) 



2k+l 



d_y n 1 

dx J cosh a; 



x=tv\\\/2U 



x=ir\\\/2U 



(Recall that throughout the paper U > 0.) Now for x > expand 1/ coshx = 
2 X]^Lo( — l) n e~( 2n+1 ) x and 2 arctan e x = 7r — 2 ^^L ^+1 e~( 2n+1 ) x . Interchanging sum- 



mations over n and and using the power series expansion of I^(z) : we obtain eq. (2.4). 

Since I v (z) ~ (2nz)~ 1 ^ 2 e z for large \z\, we see that the expansions (|2.4[) abso- 
lutely converge for |A| > 1, and in fact they are convergent also for |A|=1. We note 
that a complementary small-(A/L r ) expansion can be obtained from (|A.1|) by expanding 
1/ coshx = Xl^Lo ^tt" ™' where i? n are Euler's numbers. This way we arrive at 



2A- 



E. 



n+2k 



— 2-"--+ 1 A-!(A- t 1)! (2^) 



2k+l 



(A.2) 



k=0 



(\\ - \ ^ 1 /-7rA\ n E n+2 k-i f k 

PcW - 2^ n \ \ 2 U ^ 2 2fc (£;!) 2 \2U 
«=i v 7 fc=0 v ' 

n odd 

n = > 
tl even 

This expansion converges for |A| < U. 



In preparation to the analysis of the massive dispersion relation (|2.1j ) we need several 
definitions and lemmas. 

Definition: For m = 0, 1, 2, . . . and integer £, define the constants Ag and as the 

coefficients appearing in the expansions 



-I) 



/ 7 \ 2(m+l) 

( £ ) (v^ + T - a) = £(-l) m+ '(4m + 1 " 2£ ) !! + 1) 



(A.3) 
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where, as usual, the double- factorial stands for nil = 1 • 3 • . . . • n when n is a positive odd 
integer and (—1)!! = 1. 

The coefficients Ag 71 ^, B^ satisfy the following recursion relations (for m = 1, 2, . . .) 
and initial conditions: 

A (m) = A (m-i) + (4m _ 2t)A { ™~ X) , = B { ™~ X) + (4m + 2 - 2£)B < £ ™[ 1) ^ 

A^p = B^ = 5^0 j 

from which it is easy to see that = B^ = for £ {0, 1, . . . m}, so that the sums 

in ( |A.3| ) are in fact finite. 
Lemma 1: For n = 0, 1, 2, . . . 

^ J ' T(2m + i - £) (2n-l)!! 

t(-lY(4m + 1 - 2.)!! < )^r(2m + n + |-£) = (2m + 2n - l)!!(2m + 2n + 1)!! 
^ J £ r(2m + f-£) (2n-l)!! 

(A.5) 

Proof: Equate powers of a 2 on both sides of ( |A.3| ) after expanding them using the binomial 
expansion (a 2 + l) s = Y^Lo (n) ^ 1 ^ < ' 
Corollary: For fc G (-|,f) 



oo . 2m+l ; m 



(m) 



00 „:„2m /„ 

COS ' 



j: -20" ^ m "" 



(A.6) 



^ (2m)! 



Proof: For the first line use k = arcsin(sin/c) = J2m=o ^(2m+i)^ sin 2m+1 fc and then 
replace the numerator of the coefficient here by the lhs of the first line of (|A.5|) , with n—0. 
For the second line expand cos k = yl — sin 2 k = 1 — J2m=i ^ 2m ~ 3 |^ 2m ~ 1 ^ !! sin 2m k and 
then use the second line of ( |A.5| ) with n=0 and m replaced by (m — 1). ■ 



Lemma 2: For m = 0, 1, 2, . . ., let 



Sm{z) = Y,{- l ) lA lZ 2m -* K 2m-l{z) , T m {z) = ^(-l)^fV™ +1 -^ 2m+1 _,(z) , 
1=0 1=0 

(A.7) 



where the K u (z) are modified Bessel functions. Then 



S m (z) = z 2m K (z 



T m (z) = z^K^z) . 



(A.t 



Proof: Using ( |A.4|) and the recursion relation 8.486(10) of |?2| for the K v (z), one obtains 
S m+1 (z) = z 2 S m (z) and T m+1 (z) = z 2 T m (z) with S (z) = K (z) and T (z) = zK x (z), 
from which (|A.8| ) immediately follows. I 



Equipped with the above results, consider (|2.1| ). Expanding the sin(u;sin/c) and 
cos(ci; sin k) in power series in u sin k and integrating term by term using 6.621(4) of |2!| we 
find 



Ps (k) = k-2J2 



m=0 



{2m + 1)1 



sin 



^(-l)-i 



n=l 



, \ 2m -, 

a \ 1 



e s (k) = U - cos A; + 2 ^ 



sin- 



da J V" 2 + 1 

2m 



a=2nU 



- (2m)! 

m=0 



n=l 



da 



=2nt/ 

(A.9) 



Now using (|A.3| ) and ( |A.6|) yields, after reordering the summations over I and 



n, 



oo . 2m+l j, 771 
m=0 V £=0 



(m) 



l + 2^(-l)-(a 2 + l)- 



•(2m+|-£) 



n=l 



a=2nU 



e s {k) = [/-l + 2j](-ir 1 (v / ^TT 



Q: 



n=l 



=2nC7 



+ £ 



00 • 2m ; m — 1 

5E__* ^(-l/(4m-3-2£)!!sJ m - 1) 



' (2m)! 

00 

l + 2 ^(-l)-(a 2 + l)-( 2 — i-^) 



n=l 



a=2nU 



(A.10) 



The sums over n encountered in the last step have been analyzed in detail by Fisher 
and Barber in ||23|| . where they are referred to as "remnant functions" (they are related to 
the so-called Epstein-Hurwitz zeta function). Using the notation introduced in eq. (2.25) 
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of ||23||, we have 



n=l 



a=n/y 



3 5 
2' 2'~2 ! 



E 



Va 2 + 1 



n=l 



_y 

2 

l 

4y . 



R[,\ (y*) + 21n2 



"1/2 



(A.ll) 



The crucial step now is to use eqs. (2.26), (6.20), and (6.21) of p3|, which leads to 



ps(k) = y 



e s (k) 



oo . 2m+l ; m 

^ £ y (-l)'(4m - 1 - 2£)U A 

(2m + 1 )! ^ v ; v ; 



M 

i 



x 



2 A sin 2m+1 k 



2U' 1 



(4m - 1 - 2£)!! 



n=0 



U S (2m + 1)! E 5 ™^ 

m=0 v y n=0 



^ v ^ Ki(U n ) \ - 



oo . Om t m — 1 

sin* fc 



n=0 



^ (2m)! 



^(-l) £ (4m-3-2f)!! S 



(m-l) 



£=0 



oo 

77 Z-^ 



#l(£/r. 



(4m-3-2£)!! 



00 „• 2m t, oo 



n=0 



n=0 



U n U ^ (2m)! 



n=0 



(A-12) 



where ?7 n = (n + 5)^7, and in the second lines in each formula we employ the definitions 
( |A.7| ) . Now invoking ( |A.8|) and reordering the summations over n and m, we finally obtain 



(U- 
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